I 



in 
o 



- 1—1 

X 

5-H 



I. INTRODUCTION 



Exact transmission moments in one-dimensional weak localization and 

single-parameter scaling 

Jean HeinrichsQ 
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We obtain for the first time the expressions for the mean and the variance of the transmission coef- 
C^*) , ficient for an Anderson chain in the weak localization regime, using exact expansions of the complex 

transmission- and reflection coefficients to fourth order in the weakly disordered site energies. These 
results confirm the validity of single-parameter scaling theory in a domain where the higher trans- 
mission cumulants may be neglected. We compare our results with earlier results for transmission 
cumulants in the weak localization domain based on the phase randomization hypothesis. 
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m ' 

The validity of the single-parameter scaling hypothesis in the scaling theory of localization of Abrahams et al. Q 
has been the object of numerous discussions during the decade following its publication in 1979. This scaling theory 
assumes that the scaling of a typical conductance g as a function of size L of a d-dimensional disordered sample is a 
function of only one parameter, namely the typical conductance itself. 

Soon after the appearance of Ref.|l| it became clear that one should consider the scaling of the full probability 
distribution of conductance (or the distribution of a related scaling variable) rather than simply that of a typical 
conductance . Single-parameter scaling then requires that the distribution of the chosen variable should scale to a 
universal form depending only on one parameter. A simple example of such a statistical distribution is the familiar 
normal distribution in which the large scale dispersion a 2 decays to zero via the law of large numbers relating it to 
the growth rate of the mean value, the only parameter of this distribution. 

Many aspects of scaling of distributions of transport related quantities such as the resistance p = ~ have been 
discussed some years ago by Shapiro and coworkers Q and in a more recent review Q. 

Until recently analytical discussions of single-parameter scaling have invariably relied on the so-called phase ran- 
domization hypothesis 0- This hypothesis assumes the existence of a microscopic length scale, shorter than the 
localization length, over which the phases of the complex reflection- and transmission coefficients of a disordered 
sample are randomized in such a way that their distribution becomes uniform. 

Let us briefly recall the scaling results for a ID-disordered chain of length L. In this case a convenient scaling 
variable related to the conductance is the finite scale Lyapunov exponent 



1 , ( 1 
^ = 2L ln (! + ~g 



i3 1 l 

a; =-w 1 ^ > (!) 

o 

where the second line follows from using the Landauer formula, p = - = ^ , which relates the resistance to the 
transmission (i^) and reflection (r^) coefficients of the chain. The advantage of the scaling variable 71, is that under 
the assumption of complete phase randomization it obeys a normal distribution 0, Q of the type alluded to above, 
in the strong localization regime, L >> L c (where L~ l = 7 = limL^oo 7l is the inverse localization length), for weak 
' disorder. The distribution is centered at the asymptotic value L~ l , with a dispersion 



a 2 -var 7L --i , (2) 

i.e. it describes single-parameter scaling in terms of 7 = L~ x . The Eq. ©coincides with the general large number 
law property for the validity of single-parameter scaling (SPS) for ji proposed recently by Deych et al. for the 
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general case where phase randomization is not hypothesized. In their work |{| Deych et al. studied the condition (j2J 
analytically for an Anderson chain with a Cauchy distribution of site energies. This has led them to identify a new 
characteristic length beyond which single-parameter scaling holds true in their system On the other hand, under 
the random phase assumption one-parameter scaling also exists in the weak localization (quasi metallic) regime for 
weak disorder [3,0], where L << L c . In this case the resistance p « l(p ~ 2L7) follows an exponential (Poisson) 
distribution 0,3 depending on the sole parameter L~ x . The condition for one-parameter scaling analogous to (J2J) 
then reads (see also Sect. Ill) 



° 2 = {lL? , (3) 

where (. . .) denotes averaging over the disorder. The Eq. J3J implies one-parameter scaling in a domain where the 
higher cumulants of the resistance distribution are subdominant. We also note that is exactly verified (with 
(Tl) 2 = 7 2 ) by the exponential distribution of resistance obtained in the framework of the random phase assumption 
for L « L c . 

The purpose of this paper is to analyse the validity of one-parameter scaling theory in the weak localization regime, 
by deriving exact analytical results for the mean value and variance of the transmission coefficient th for a weakly 
disordered Anderson chain, without using phase randomization. The resulting expressions for the variance a 2 and 
the mean value (tx) of the variable are shwon to obey approximately (for L » a, with a the lattice parameter) 
a universal relation analogous to ©. Here the mean (tl), which turns out to depend on L is completely defined 
in terms of a fixed localization length. The complex transmission coefficient T^it^ = \Ti\ 2 and reflection coefficient 
Rl(tl = \Rl\ 2 ) are obtained in Sect. II by iterating the exact coupled recursion relations obeyed by these quantities 
[a,l3i f° r weak disorder in the regime L << L c . These results are used for obtaining explicit forms of the first two 
moments of the transmission coefficient tt and of the corresponding Lyapunov exponent 7/,. These moments are 
discussed in connection with the validity of one-parameter scaling theory in Sect. Ill, where they are also compared 
with previous results based on the random phase assumption. 



II. TRANSMISSION IN THE WEAK LOCALIZATION REGIME 



Consider a disordered linear chain of length L = Na composed of N disordered one-orbital sites of spacing a = 1 
described by the Schrodinger equation 



(fn+l + (fin-1 + £n<-Pn = Elf n , (4) 

where ip m is the wavefunction amplitude and e m the random energy of the site m measured (like E also) in units of the 
constant hopping rate between nearest-neighbour sites. We assume the e m to be identically distributed independent 
gaussian variables with mean zero and correlation 



(^m £n) — ^Q^m.n ■ (5) 

The usual Anderson model of uniformly distributed site energies between — W and W correspond to Eq = The 
disordered chain is connected to electron reservoirs by semi-infinite non-disordered leads defined by sites m > N and 
m < 1, respectively. As indicated in Sect. I, we whish to calculate the transmission coefficient = \Tn\ 2 and to 
study its statistical moments in the weak localization regime L « L Cl for weak disorder. The inverse localization 
length is defined by 



7 = — = — lim lniiv . (6) 

The complex transmission- and reflection amplitudes for a plane wave incident from the right (n > N) with 
wavenumber-/c and energy 



E = 2cosk, < k < 7T 



(7) 
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obey the exact coupled recursion relations 0, 



e ik T N . 
l-iv N (l + e 2ik R N -i) 



Tn - i 77. n , o2 ifcp \ ' ( 8 ) 



e 2ik R N -! + WnQ- + e^ iW) 
ii/iv(l + e 2lk R N - 



^ = ^T . ( 10 ) 
2smfc 

which relate the amplitudes Tn(Rn) for the chain of N disordered sites to those for a chain of N — 1 sites, with 
the initial conditions To = 1 and Rq = 0. The choice of wavenumbers in (JJJ is such that the leads eigenfunctions 
tf n ~ e ±lfc " correspond to Bloch waves traveling from left to right and from right to left, respectively. The solution 
of jnj) may be readily obtained in terms of the reflection amplitudes for chains whose lengths are comprised between 
1 and N - 1: 



n=l 



For weak disorder we expand to successive orders in the site energies e n . As will be readily seen below the 
validity of this perturbation expansion is restricted typically to chain lengths L << L c (the weak localization regime), 
where L c is the localization length for weak disorder first derived by Thouless @ . The study of the dispersion of the 
transmission coefficient requires to include terms up to fourth order in the expansion of We thus obtain 



N N 



T N = e* kN [1 + ]T (P m - 1) + J2 ( P « - !)( P « - !) 
m=1 to, n = 1 

m ^ n 

+ (Pm-l)(P»-l)(P p -l)+ ]T (P m -l)(P w -l)(P p -l)(P 9 -l) , 

m,n,p m^n^p^q 
m ^ n ^ p ^ m m ^ n ^ p ^ q 



n 7^ q 7^ m ^ p 



1 



P j = — — , Xj - -i Vj (l + e^Rj-i) . (12) 
i -t- 



In each one of the successive terms in (1 1 21) , we then retain contributions up to fourth order in the variables Vj , via 
the expansion of the Pj in power of Xj and corresponding expansions of Rj in successive order contributions in the 
reflection amplitudes, 



Rn=Y. R n ] ■ ( 13 ) 
1=1 

The contributions R^ at successive orders q are determined by the perturbation equations obtained by expanding 
the recursion equation analogous to 10 relating R n and R n -i for systems of length 1 < n < N — 1 and n — 1: 

(14a) 

'I , (14b) 

',, „4»&d(1) 2 n,2 2ifcp(l) -3 
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whose solutions are 



p) 



P =i 



^ 2) = -E^ 2lfe 



(n-p) 



a 2ik(jp-p) 



q=l 



P =i 



(15a) 
(15b) 
(15c) 



Next, by using H12I13JI and (15.a-15.c), we express the transmission coefficient — TV | 2 to fourth order in the 
site energies and, finally, we find its average (t^) using together with the factorization property for independent 
gaussian variables, 



(^m^n^-p^-q) — {^m^n) (^pE-q) ~t~ (^-m^-p) (&Ti£q) i^-m^q) (^-n^p) 

which is valid for arbitrary indices m, n,p, q different or not |9j. 

These straightforward but rather tedious calculations yield the final result 



(16) 



(t N ) = 1 - 



2N 2N 2 



1 



L c LI 

which is exact to quadratic order in e\. Here 



o — — - 



2 sin 2 kN / cos 2 kN 
cos^ k 



N N 2 sin 2 k 



+ 0{et) 



(17) 



1 



p 2 
-0 



F 2 



\sin 2 k 2(4 -E 2 ) 



(18) 



is Thouless expression |8j for the localization length, which is valid to lowest order in the weak disorder. The 
perturbation expression Ijl7|l shows that besides the smallness of Eq its validity requires typically N << L c , which 
implies that the eigenstates are delocalized on the scale of the chain length. We further note that while the lowest 
correction in (t^) for weak disorder is of order the next higher contribution proportional to £q includes three 



types of terms respectively proportional to (77-) , ("T") £ o an< ^ £ o- By a similar calculation of the second moment, 



(t%), we obtain 



(t 2 N ) 



47V 4iV 2 



1 sin 2 fc7V ( 2 cos 2 kN sin 2 kN 
4 



N N 2 sin 2 fc 



cos 2 k 



2 sin^ k 



(19) 



The calculation of (t 2 N ) has been greatly facilitated by making use of the expression of tN to quartic order and of the 
corresponding average l|17f) . From l|19[) and (|18f) we finally get 



var tN = 



12N 2 



2-1- 



1 sin 2 kN ( . sin 2 kN 



N 3N 2 S in 2 fc 



2 sin 2 k 



+ O(e 6 ) . 



(20) 



The first and second moments and the variance of in are easily found from l|17f) and l|19[) by expanding In t n 
in powers of the small quantity 1 — ti\r- Thus we obtain 



(In) 



2N 



(lnt N ) 



1 



AN 
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N 



1 sin 2 kN 



2N 2 



?k 



cos 2 kN sin 2 kN 



cos 2 k 



4sin2fc 



o(4) 



(21) 



(iD = J2 



3 1 sin 2 fc7V 
~ N ~ N 2 sin 2 k 



sin 2 kN 
2 sin 2k 



■o(4) 



(22) 
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a 1 = var -f N = — 



1 1 sin 2 kN f sin2fciV\ 
~N~7P sin 2 k V + 2sin2£; ) 



(23) 

By ignoring the terms of orders -h and relative to unity in the square bracketts of (|21I23|I it thus follows that 

* 2 = 6((JN)) 2 + O(e 6 ) . (24) 

III. DISCUSSION AND CONCLUDING REMARKS 

The equation (|24fl is a universal relation expressing the variance of the scaling variable in terms of the mean, 
which is itself defined in terms of the Thouless localization length (|18fl alone, when terms of order j? and are 
neglected relative to unity in (|21(l . It is natural to expect that, within the same approximation, similar universal forms 
in terms of (7_/v) exist for the higher fluctuation cumulants of jn. Thus our work suggests that the exact distribution 
of 7jv (which does not rely on phase randomization) in the weak localization regime obeys single-parameter scaling. 
Note that l|24|l differs by a factor of 6 from the universal relation © for the weak localization regime obtained under 
the phase randomization assumption. We observe in passing that the mean of 7jv in Ij21|l differs from the Thouless 
inverse localization lengths (|18|l . This is because the inverse localization length is by definition, the non-random 
central limit of jn. In order to obtain L^ 1 from our analysis for TV << L c one has to let |eo| — > before taking the 
large N limit. In this limit only the leading order term of the weak disorder expansion of (tjy) in l|17fl is relevant 
for finding L~ l . We believe that the reason why this term yields the correct (Thouless) localization length in spite 
of the fact 7 at is not self-averaging in the domain N « L c is because the non-self averaging is only marginal i.e. 

2 

limjv^oo hm| eo |_>o j^jj ^constant in this case. 

The aim of the present work has been to discuss the scaling of -yj, in the weak localization regime without using the 
phase randomization hypothesis. In this context we now compare the above results for the transmission coefficient 
moments with earlier results obtained from an invariant imbedding model [H 13 analysed by assuming phase 
rendomization llOf. The coupled invariant imbedding equations for the complex reflection (Rn)~ and transmission 
coefficients (TwYpllL IT^ have since been shown |fj to correspond to the continuum limit of the Schrodinger equation 
for the one-dimensional Anderson model for weak disorder. This makes the comparison of the present results with 
those of ref.^3 the more relevant. Under the random phase assumption we obtained the following exact results 0] 
for the transmission coefficient moments for L << L c 

iW- 1 



W)=[l + 2m T , m = 0,1,2,3,... , (25) 
which leads to the distribution of the transmission coefficient 01 

P ^ n) = ^ (W2A,_1) > ( 26 ) 

which depends on the single scale parameter L c . Note that the exact moments (|17|) and i|19|) which do not rely on 
the random phase assumption reveal important deviations from Eq (|25|l at order even when the terms in and 
in the square bracketts arc neglected 14]. The exact random phase expression for var t^, 

47V 2 

va,rt N = — , (27) 

obtained from (|25() differs also from the exact expression l|20l) to order jj. For completeness sake we finally list the 

expressions for the first and second moments and the variance of obtained from expanding the logarithm of 125|) 
for m = 1 and m = 2 through order Eg: 
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<7^> 
° 2 = ({In)) 



LI ' 
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(29) 
(30) 
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